o 
o 

(N 
> 

o 

Continuity corrections for certain perpetual Annerican 



> 

o 



in 
o 

B 

> 

X 



and Bermudan options on multiple assets 

Frederik Herzberg'^ ^ 



^Keywords: American options, barrier options, exercise regions, continuity 
corrections. 

^AMS Mathematical Subject Classification (2000) : 60G51; 60G40; 91B28. 
^Abteilung fiir Stochastik, Institut fiir Angewandte Mathematik, Universitat 
Bonn, D-53115 Bonn, Germany (herzberg@wiener.iam.uni-bonn.de) 

''Mathematical Institute, University of Oxford, Oxford 0X1 3LB, England 



Abstract 



In a general Feller martingale market with several assets, the existence of optimal 
exercise regions for multi-dimensional Bermudan options can be established by 
reference to Neveu's theory of Snell envelopes - and also, as will be shown, more 
directly from standard martingale arguments combined with the strong Markov 
property. Based on this, in the framework of a log-Levy martingale market, 
explicit formulae and asymptotic results on the perpetual American-Bermudan 
(barrier-like) put option price difference ("continuity correction") near the exercise 
boundary will be proven, under the - of course, fairly restrictive - assumption 
that the logarithmic optimal exercise region, subject to the barrier, does not 
depend on the time mesh size and is, up to translation, a half-space. 

For this sake, Wiener-Hopf type results by Feller will be generalized to higher 
dimensions. It will be shown that an extrapolation from the exact Bermudan 
prices to the American price cannot be polynomial in the exercise mesh size in 
the setting of many common market models, and more specific bounds on the 
natural scaling exponent of the non-polynomial extrapolation for a number of 
(both one- and multi-dimensional) market models will be deduced. 
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0.1 Introduction 



American options are financial securities that can be exercised at any future 
tinne before maturity, or, in the case of perpetual options, at any date in the 
future. Bermudan options, on the other hand, can only be exercised at certain 
dates in the future (and therefore are sometinnes also referred to as discretely 
sampled American options). An option on a multiple asset {or: multi-dimensional 
option) is a contract whose payoff function depends on more than one asset. An 
(American or Bermudan) barrier option can only be exercised if some asset price 
falls below or rises above a certain level. 

What the holder of an American or Bermudan option can gain from exercising 
the option depends on the exercise strategy she chooses. The issuer of the option 
has to set as price the least upper bound of all expected gains from such exercise 
strategies. 

In a memory-less market, this supremum coincides with the supremum over 
all exercise strategies that only look at the current price, and this supremum, of 
course, is a maximum, since at each exercise time one will only have to check 
whether the payoff one would get from exercising now is still less than the current 
option price (which in turn is taken to be the supremum mentioned above). By 
choosing smaller and smaller steps between the (discrete) exercise times of a 
Bermudan option, one can approximate an American price. 

Whilst the reasoning of the previous paragraph has long been brought to 
mathematical precision for one-dimensional options, and nowadays seems to be 
fairly well-understood for higher dimensions as well, it appears still to be quite 
difficult to find references that rigorously prove the existence of exercise regions 
for multi-dimensional Americans or Bermudans. Therefore, this paper fills a 
gap in the existing literature by providing a rigorous introduction to (multi- 
dimensional) American and Bermudan option pricing in Section U 

We will thus first of all give a rigorous account of the existence of optimal ex- 
ercise regions for Bermudan options on multiple assets, their coinciding with im- 
mediate exercise regions and (as an immediate Corollary) their time-stationarity 
in the case of perpetual options. Explicit characterisations of immediate exer- 
cise regions for (non-pereptual) Bermudan options can be found in the work of 
Broadie and Detemple [8j, Ekstrom [I4j as well as Peskir [35], whereras the 
immediate exercise boundary of one-dimensional perpetual Bermudan options is 
discussed in an article by Boyarchenko and Levendorskii (In the spirit of 
Boyarchenko and Levendorskii's paper, we shall look at the operator equations 
corresponding to - multi-dimensional - perpetual Bermudan pricing problems in 
an Appendix.) 

Having proven the existence and time-stationarity of the exercise regions for 
Bermudan options, we shall then move on to studying continuity corrections - 
that is, the differences of an American price and a Bermudan price of a given 
exercise mesh size - of a perpetual put option. 

In particular, we shall derive asymptotic bounds on these corrections (con- 
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ceived of as a function of the Bermudan's exercise mesh size) as the Bermudan 
exercise mesh size tends to zero. For this purpose, a Wiener-Hopf type result of 
Feller's [TF, p. 606, Lemma 3] will be generalised to higher dimensions. 

The motivation for finding such continuity corrections is, of course, to be 
able to extrapolate from a finite number of Bermudan option prices to an approx- 
imation to the American price (both with the same particular payoff function). 
Therefore, the problem has attracted interest for quite some time, and computa- 
tionally useful results for one-dimensional options have been established as well. 
One can find these scaling results in the works by Broadie, Glasserman and Kou 
[To] , by Howison (23 as well as by Howison and Steinberg [22j. 

The class of payoff functions g to which our continuity correction results 
will be applicable is the class of half-space barrier-like payoff functions - given a 
discount rate r > 0, the barrier-like payoff functions are defined as exactly those 
payoff functions with the property that for sufficiently small excercise mesh sizes 
s < So. the optimal logarithmic exercise region G := Gg does not depend on s. 
If this G is, up to translation, a set that is closed under addition in with its 
complement also being closed under addition, then the payoff function will be 
called a half-space barrier-like payoff function. We will also show that all convex, 
closed and +-closed sets with on the boundary whose complements are also 
+-closed, are merely half spaces in the following sense: If i7 C M"' is closed, 
convex and +-closed with CH also being +-closed, and G dH, then there is a 
VH e M'^ such that H = {xeR'^ : *x • y// > O}. 

Note that the price of a (possibly multidimensional) knock-in barrier option 
with payoff function g where the possible exercise region G determined by the 
barrier(s) is a subset of the optimal exercise region of the corresponding Amer- 
ican option will be the same as the price of the American/Bermudan option 
with barrier-like payoff function gxG, since the option will be exercised immedi- 
ately after the barrier has been hit. This vindicates the term "barrier-like payoff 
functions". 

This paper is based on Chapter 1 and Appendix A of the author's thesis [19] . 
0.2 Notation 

We are following largely standard probabilistic notation, as can be found for 
instance in the works by Ito and McKean jr [27] or Revuz and Yor [37]. In 
particular, if P is a probability measure on a cr-algebra C, Z a random variable 
and C G C, thenE[Z,C] = f^ZclF, and ifP[C] > 0, thenE[Z|C7] = j,E[Z,C], 
whilst the conditional expectation of Z with respect to a sub-a-algebra Cq C C 
will be denoted by E [Z\Co] 

Both A O B and A C -B for sets A and B will mean that A is a subset of B 
(possibly A = B), whereas we shall write A B to express that ^ is a proper 
subset of B. 

Finally, for any subset A C M'', shall denote its complement := 
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1 Definitions and basic facts on Bermudan and 
American options 

1.1 Terminology 

Our first definition is a notational convention. 

Definition 1.1. Let d G N. By exp : M>o'^ and In : R>o'^ R'^ we 

denote componentwise exponentiation and taking natural logarithms compo- 
nentwise, respectively. 

Remark 1.1. For any d ^ 'N, is a Lie group with respect to componen- 
twise multiplication ■ : {x,y) (2;jyi)ie{i....,d}- algebra is the vector 
space with its usual (componentwise) addition. The exponential map from 
the Lie algebra (M'^,+) into the Lie group (M'^,-) is componentwise exponen- 
tiation exp : X I— > (e^')-g|^ Therefore the abbreviation introduced in 
Deftnition \Ll\ is consistent with standard notation. 

Definition 1.2. Let T be a positive real number. Consider a real-valued 
stochastic process X := (-^t)tg[o,T]; adapted to a filtered probability space 
(rj, (^()tg[o,T]; -P) • will call X a logarithmic price process for an asset 
with continuous dividend yield 5 (for short, a logarithmic price process or 
simply log-price process^, if and only if there exists a probability measure 
Q equivalent to P on J-t and a constant r > such that the stochastic 
process exp {Xt — rt + (5i)jg[oT] ^•^ ^ martingale with respect to the filtration 
T := {^t)t<=lQ,T] o.iT'd the probability measure Q. In this case, such a Q is 
called a martingale measure and r a market price of risk or a dicsount rate 
for the stochastic process X and the probability measure P. 

Definition 1.3. Let d £ N. A d-dimensional basket is a d-tuple of loga- 
rithmic price processes such that there exists a probability measure Q and 
a market price of risk r > such that Q is a martingale measure and r a 
market price of risk for all components of the d-tuple. 

For the rest of this section, we will adopt the ternninology and the notation 
for Markov processes of Revuz and Yor |37] . 

In particular, for all probability measures v on B{W^), Qy is the probabil- 
ity measure induced by the transition function {Qs)s>o via the lonescu-Tulcea- 
Kolmogorov projective limit construction, cf. Revuz and Yor Theorem 1.5]). 

For any d G N, we will denote the cr-algebra of Borel subsets of by 
B iW^) . 

Definition 1.4. Let again d G N. family Y := iY.'')^^^^ of W^-valued 
homogeneous Markov processes Y^ adapted to a filtered probability space 
(^,J-^Q^ with respect to T, with transition function (Ps)^>q and initial 
measure 6x, is called a d-dimensional Markov basket with dividend yields 
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6i,...,5d > if and only if there is a homogeneous transition function 
{Qs)s>o measurable space (M.'^ , B (M.'^)'j and a constant r > such 

that the following three assertions hold: 

1. The process is a Markov process with transition function (Qs)s>o 
with respect to T for all x G W^. 

2. The process i{ex.p{{Y^t)- — {r — 6i)t))-^^^ is a martingale 
with respect to J- and Qs,^. ■ 

3. The measures Ps^ and := Qs^ are equivalent for all x S M'^. 

In this case, {^^)x^^d is called a family of martingale (or: risk-neutral) 
measures associated with Y, and r is called the discount rate for Y. 

The expectation operator for the probability measure will be denoted 
by for all x G M^. 

// the transition function P is a Feller semigroup, then we shall refer to 
Y as a Feller basket. 

If P is a translation-invariant Feller semigroup, we shall call Y a Levy 
basket . 

Remark 1.2. A priori, it is not clear if there are logical connections between 
the three assertions in the previous Definition \1.4[ in particular the author 
does not know whether the third assertion implies the first one. 

Notational convention 1.1. // no ambiguity can arise, we will drop the 
superscript of a Markov basket. Thus, in the notation of Definition \l.^ we 
set 

E-[/(y.„...,y.j|^,] :=E^[/(y,^,...,y,:)|^,] 

for all s > 0, n Gf^ and n-tuples of stopping times r = (ti, . . . , r„) whenever 
/ : M" ^ M is nonnegative or f (y^^, . . . , Y^^J G (P^). Here we are using 
the term "stopping time" as a synonym for M.-^- -valued stopping time, that is 
a stopping time with values in [0, +oo] . 

Also, since we are explicitly allowing stopping times (with respect to the 
filtration generated by a process X) to attain the value +oo, we stipulate 
that the random variable f (Xr) (for any Lebesgue-Borel measurable function 
f) should be understood to be multiplied by the characteristic function of 
the event {r < +00} . Formally, this can be done by introducing a constant 
A W^, called cemetery, and stipulating that Xr = A on {t = -\-oo} and 
/(A) = for all measurable functions f (cf. e.g. Revuz and Yor pp 
84,102]). 

We will not formally define what we nnean by an option itself, but rather 
provide definitions for the concepts of expected payoffs and prices for certain 
classes of financial derivatives. 
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Definition 1.5. Consider a d- dimensional Markov basket Y with an asso- 
ciated family ¥' of martingale measures and discount rate r > 0. 

The expected payoff of a Bermudan option with (log-price) payoff function 
g : R'^ ^ M>0 on the underlying Markov basket Y with exercise times in 
J C [0, +oo), log start-price x and maturity T G [0, +oo] is defined to be 



U'{T){x) := Ug\T)ix) := sup E^' e-^^^^'^ 'g (X.^t) 

T stopping time, r(n)CJu{+oo} 

The expected payoff of a perpetual Bermudan option is the expected payoff 
of a Bermudan option of maturity +oo. 

The expected payoff of a Bermudan option with exercise mesh size h > 
is the expected payoff of a Bermudan option with exercise times in h -Nq . 

The expected payoff of an American option is the expected payoff of a 
Bermudan option with exercise times in [0, +oo). 

We shall call the expected payoff of a Bermudan option (or an American 
option) a Bermudan option price (or an American option pricej if and only 
if the martingale measures associated with the underlying basket are unique 
(that is, if the market model described by P, T and X is complete^. 

In recent years, there has been increasing interest in incomplete nnarket nnod- 
els that are governed by general Levy processes as log-price processes, as is not 
only witnessed by a tendency in research papers to focus on Levy process settings 
(for instance Boyarchenko and Levendorskii [7j; Asmussen, Avram and Pistorius 
f4j; 0ksendal and Proske [33J, to take a random sample). Even textbooks, 
such as Karatzas' [28j and Mel'nikov's [31] introductory works, are putting con- 
siderable emphasis on incomplete markets. Finally, "Levy finance" has already 
been treated in survey articles intended for a general mathematical audience, 
e.g. Applebaum's article [3J. 

Whilst quite a few of our results will apply only to the Black-Scholes model, 
some of our arguments also work for market models where the logarithmic price 
process merely needs to be a Levy process. However, the theory of incomplete 
markets (and all Levy models other than the Black-Scholes model are incomplete) 
is not the focus of this thesis, therefore we did not endeavour to go beyond the 
Black-Scholes model where this caused technical difFiculties rather than making 
proofs easier to read. 

Whilst there are some points to be made about market failures on stock 
markets that might entail arbitrage opportunities (for example, when assets are 
traded simultaneously on several stock exchanges, or in the event of insider 
trading), the transaction costs to exploit these arbitrage opportunities usually 
tend to be close to the actual gain that can be achieved through taking advantage 
of the arbitrage. Therefore we shall, for the sake of mathematical simplicity, 
merely refer to the works of Corcos et al. [12|, Imkeller [23 as well as Imkeller, 
Pontier and Weisz [S^. and impose a strict no-arbitrage assumption - which 
under certain regularity conditions on the basket is equivalent to the existence 



,-r{T^T) 
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of an equivalent martingale measure (a measure equivalent to the market model 
under which the discounted asset prices are martingales), cf. Karatzas [2H1 
Theorem 0.2.4] and references therein. 

Example 1.1 (A few common examples). 1. The price of a Euro- 
pean call option on a single asset with maturity T and strike price K 
is the price of a Bermudan option with the set of exercise times being 
the singleton {T} and the (log-price) payoff function (exp(-) — K) V 0. 

2. The price of a perpetual American put of exercise mesh size h > on 
the arithmetic average of two assets in an underlying basket with strike 
price K is the price of a Bermudan put option with the set of exercise 
times being the whole of the half-line [0,+oo), the maturity being T 
and the payoff function [k - ^M{-V)+^^M{-)2) S^ y g. 

3. Consider a perpetual Bermudan call option on a single asset that con- 
tinuously pays dividends at a rate 5 and whose logarithm follows a 
Markov process Z adapted to some probability space {Q, {J^t)t>o, P) ■ 
Then, in order to exclude arbitrage, we will have to require the ex- 
istence of a family of measures F' such that each is equivalent to 

(in particular, Pf^^ = 5x) and such that (e"''*'*'''*'*'^') is a F^- 

martingale for all x G M°'. The expected payoff of the option will then 
be 

[/'^•^"(•) = sup E- [e-"" (e^- - K) V O] 

T stopping time, r(f2)C/iNoU{+oo} 

1.2 Convergence of Bermudan to American prices 

In this paragraph, we shall give a formal proof in a general setting that the 
price of a Bermudan option with equidistant exercise times (of mesh h) before 
maturity indeed converges to the American price as h tends to 0. 

As an auxiliary Lemma, let us remark the following elementary observation: 

Lemma 1.1 (Lower semi-continuity of sup). If I is a set and 

{'^k,t)i^i fcgf^g is a family of real numbers, then 

supliminf Ofc^^ < liminf supa^^^. 
Proof. We have the trivial estimate 

for all A: G No and £o G I , therefore for all n € N and £o S I, 

inf supafc,^ > inf ak^i^, 

k>n n k>n 
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thus 



hence 



inf supak,e > sup inf ak^, 

k>n f Pn k>n 



sup inf supak/ > sup sup inf ak/^ = sup sup inf afc/g. 

This is the assertion. □ 

This estimate enables us to prove the following Lemma that is asserting the 
approximability of expected payoffs or prices of American options by sequences 
of expected payoffs or prices of Bermudan options, respectively. 

Lemma 1.2. Let d £ N, T > (T = +00 also possible), x G W^, and let 
X be a d- dimensional Markov basket. Suppose the payoff function g > is 
bounded and continuous, and assume X has a modification with continuous 
paths. 

Then, if the American expected payoff U^^'~^'^\T)(x) is finite, one has 

u^^'+°°\t)(x) = nmC/™o(r)(x) 

hlO 

= supU^''^''{T){x). 
feeN 

Proof. Consider a sequence {hk)ke.no £ (I^>o)^" such that /i/j J, as A; ^ 00. 
Choose a sequence of stopping times (r^j^gp^g such that for all x G M'', 



sup 

r stopping time 



and define 

r^,fe := inf {t G /ifcNo : t > t^} . 

Then, due to the continuity conditions we have imposed on g and on the 
paths of (a modification of) the basket X, we get 

supe-^(-^^^)5(^r,AT) = sup lim e-H-^'"^^)^ (X,, ^^.t) 
and hence by the lower semi-continuity of sup, one obtains 

supe-'-(-^^^)5(^r,Ar) < liminf supe-K-^.'^AT)^ (X,,,at) ■ 

Now we can use the Montone Convergence Theorem and Lebesgue's Domi- 
nated Convergence Theorem (this is applicable because of the boundedness 
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of g) to swap limits/suprema with the expectation operator. Combining this 
with the specific choice of the sequence (T^)^gj^ , this yields for all x G W^, 



sup e^^'(^^^)(7(X^AT) =supE^ 

T stopping time '- i 



re AT) 



sup e 



< liminfE^ 

fc— +00 



< E'' 



liminf supe-K-^^''^^)^ (X,^ ^^r) 

k^oo £ ' 



liminf supE^ \e-''(^'-^^^) g (X^,,,at) 



fe— >oo 

< liminf sup E^' 

*=^°° T{Q)C/ifcNoU{+oo} 



-r(TAT) 



< limsup sup E' 

fc-*oo r(C)C/n.NoU{+oo} 

< sup E'' 

T stopping time 



g {XrAT 

e-^(^^^^g{X^^T) 



This finally gives 



sup E^' 

T stopping time 



-r(TAT) 



Hm sup E^ 



g (XrAT) 



Since the left hand side does not depend on {hk)k, we 
conclude that Wm.hiQU^^'^ {T){x) exists and is equal to 

sup. .topping time [e'''^^^^) <7 (^tAt)] • □ 



2 Exercise regions 

An American/Bermudan option price coincides with the payoff that is expected if 
one exercises at the first possible entry of the log-price process into the immediate 
exercise region (which is a subset G C M'' in case of a perpetual option and a 
subset Q gW^ X [0, +cxd) for a non-perpetual). The immediate exercise region 
for a Markov basket with payoff function g is defined as 

F-^.^ := E M'^ X J : {T - t)(x) < g{x)^ 

for a non-perpetual option with maturity T and a set of exercise times J, and 
as 

F^:=\xeM.'^ : C/™«(+oo)(x) < ^(x)} 

for perpetual options with exercise mesh size h; its optimality for a large class of 
options was established in the theory of optimal stopping (cf. e.g. Neveu [22l 
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Proposition VI-2-8] or El Karoui [El Theoreme 2.31], whose results need to be 
applied to the corresponding space-time Markov process) using the so-called Snell 
envelope (cf. Griffeath and Snell [18j). In this Section, in addition to setting 
up notation, we shall give alternative elementary proofs for the optimality of 
immediate exercise regions for put options with a discrete set of exercise times. 
Explicit characterisations of immediate exercise regions for certain special cases 
have been proven in recent years: by Broadie and Detemple [H], Paulsen [34| . 
Ekstrom [14] as well as Peskir ^35j . 

For the rest of this Section, we would like to restrict our attention to hitting 
times where, no matter if the paths of the process are right-continuous and 
the target region closed, the infimum in the definition of a hitting time is always 
attained. We will do this by imposing the condition that the range of the stopping 
time be discrete in the following sense. 

Definition 2.1. A subset I of a topological space {X,T) is called discrete 
(with respect to T) if for all x ^ I there exists an open set U 3 x such that 
inu = {x}. 

Given any discrete subset J of a M+, it is, by density of the rationals in M, 
possible to find an embedding of T into Q-|-, thus I must be countable. By 
an analogous argument, all discrete subsets of separable metric spaces must be 
countable. 

A subset of a discrete set (with respect to a topology T) is again discrete 
with respect to the same topology T, and if T comes from a linear order, the 
infimum of any discrete set is attained and therefore by definition a minimum. 

Definition 2.2. Given a discrete subset I C [0, +oo) and a Lebesgue-Borel 
measurable set G C M°', often referred to as exercise region, we define the 
stopping time 

Tq := mm{t G T : Xt G G} , 

(the superscript will be dropped when no ambiguity can arise) which is just 
the first (nonnegative) entry time in Z into G. If Q is a subset of space- 
time, that isQ <ZW^ >^ [0, +cxd) rather than space (ie W^) itself, we use the 
space-time process rather than just the process itself to give an analogous 
definition: 

r|:=min{t : {Xt,t) ^ G} , 

where to is the time- coordinate at which the space-time process was started. 
Also, for h > we set 

h hn hN 

■— ' ■— 

to denote the first positive entry time in KHq into G or Q , respectively, whilst 
finally fg := r™" and Tq := will denotes the first nonnegative entry 
time into Q and G, respectively. 
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For convenience, we will also adopt the following convention for this section: 

Definition 2.3. Let Z C [0,+oo). A stopping time t is called I -Yalued if 
the range of t, denoted by ran t, is a subset o/XU {+00}. 

The following Lemma, as well as its Corollary can be proven easily by resorting 
to the well-understood theory of optimal stopping and Snell envelopes, cf. e.g. 
Neveu and El Karoui. Our proof will be elementary. 

Lemma 2.1. Consider a discrete subset I C [0, +00). Let X be a d- 
dimensional basket with an associated risk-neutral measure P and discount 
rate r > 0. Suppose g = {K — f)\/ and e'"^' f{X.) is a ¥-submartingale. For 
all I-valued andF-almost surely finite stopping times r there is a space-time 
region B = Uueji^} ^ such that 



E 



-r(rAT) 



< E 



for all T G [0, +00) where 

f := t| = inf {uel : X„ G B^} 

and T > F-almost surely. If the set 

{e~'^^g (Xy) : v I-valued stopping time^ is uniformly ¥-integrable, then 
the latter inequality will also hold for T = +00. 

The Lemma holds in particular for I = sNq for arbitrary s > 0. 

The condition of {e~^'^g (Xy) : v X-valued stopping time} being uniformly 
P-integrable is what is known in Neveu's terminology [32. e.g. Proposition 2.29] 
as (e-™(7 ^.^^i^^^ ^^^^^-^^ ^-^^ being of class (D). 

Proof of Lemma Wl\ Firstly, we will treat the case of T < +00. Define 

Vt G T Bt := Xr ({r = t]) C M.'^ . 

Let us first of all assume that 

G J {f = t} n {r > T} = 0. (1) 

and let us also for the moment suppose 

Vt G T g{Xt) > a.s. on {Xt G Bt} . (2) 

Both of these assumptions will be dropped at the end of the proof for the 
case T < +00 in order to show the Lemma in its full strength. Now, from 
equations Q and ^ one may derive 



E 



g{XrAT),{r = t} 



E[e-'^g{Xr),{r = t}] 

■{K-f{X^)),{f = t}] 



E [. 
E 



-r(TAT) 



iK-f{XrAT)),{r = t} 



(3) 
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for all t G [0, T] n I. 

Furthermore, observe that r > f a.s. Using Doob's Optional Stopping 
Theorem (see e.g. Varadhan jlHl Theorem 5.11]), we infer from our assump- 
tion of e~'' /(X.) being a P-submartingale with respect to the canonical filtra- 
tion the assertion that (e~™/ at tat} ^ P-submartingale with 
respect to the filtration {J^tat,^tat}- Hence, if we combine this with equa- 
tion (El and note that {f = t} = {f AT = t} e TrKT for all t G [0, T) n T, 
we obtain for every t S [0, T) n X, 



E 

= E 

< /ST-E 

< -E 

= E 



rirKT) 



<7(X,AT),{f = t} 



TirhT) 



(K-/(X,AT)),{f = t} 



-■r(rAT) 



,-r(f AT) 



E 
E 



,-r(fAr) 



/(X,AT)),{f = t} 



^(^^^)<7(^fAT)),{f = t} 



(4) 



On the other hand, since f < r, if f > T, then also r > T, entailing 

f NT = T = T NT OVL {f>r}. 
Summarising these last two remarks, one concludes 



rATj 



^E[e-^'(-^^)5(X.AT),{f = t} 

5^ E[e-'-(-^^)5(^.AT),{f = t} 
texn[o,T) 

+ E 



E 



-r(rAT) 



tGXn[T,+oo) 



< 



iGXn[0,T) 



E 



~r(f AT) 



texn[T,+oo) 
^E[e-'-(^^^)5(^fAT),{f = t} 



AT J 



(5) 



In order to complete the proof for the case of T < -|-oo, let us show that 
the assumptions ^ and Q are dispensable. 



12 



We first of all simply define the stopping time 



^' •= ^C(Utei{r=t}n{9{X0>0}) ■ °0 + X{r=t}n{g{Xt)>0} ■ t 

and based on this definition, we would set 

VtGT B't := X^, {{t = t}) . 

Then 



i?; = Xt ({r' = t}) C Xt i{g{Xt) > 0}) C {5 > 0} 



hence 



yt£l g{Xt) > a.s. on {Xt G B'^} 
However, in any case 



(6) 



E 



E 



'r(TAT) 



Now, suppose at least one of the conditions Q and Q was not satisfied 
(if 1I2II holds, one may even replace r' by r in what follows). In this situation 
we consider the stopping time 



X{r'<T}u({T'>T}n{T^,>T}) + ^C({r'<T}u({r'>T}n{T^,>T})) 



00. 



If one now defines 



ytei B'l := Xr" {{t" = t}) 



then 



Vt e J n [0, T] B'l = X^u ({r' = t" = t}) = Xr' {{t = t}) = B[ 



and 



ytein (r, +00) b[' = x^" {{t" = t}) = x^" {{t" = t' = t}n {r|, > r}) 

C Xr,{{T' = t})=Bi 



hence 



Vt E J g{Xt) > a.s. on {Xt G B'/} (7) 
(because of © and we have just seen B" C B^ for all t Furthermore, 



'"S" - ^{r'<T}u({r'>T}n{r^,>T}) ' '^B' + ^C({r'<T}u({r'>T}n{r^, >T})) ' 

= X{t"<t} ■ + X{+oo>T">T}n{ri >r} ■ ^B' 
+>^C({r'<r}u({T'>T}n{r^,>r})) " °° 



oo 
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(the first hne because of X^" (C ({r' < T} U ({r' > T} n {r|, > T}))) = 0). 
Therefore 

r|„ > T on {r" > T} , 

thus (P) holds for r^„ instead of f and r" instead of r. But we have al- 
ready proven (HJ. Therefore, analogously to the derivation of © under the 
assumptions of both ^ and (P), we get 



E 



-r{T"AT) 



g{Xr" 



AT) 



< E 



On the other hand, however. 



E 



E 



-r(T'AT) 



at; 



(as r" = r' on {r < T}, as well as r" > T on {r > T}, thus t" AT = T 
r A r on {r > r}) and we have already seen that 



E 



Finally, 



E 



e-'-(^'"^)<7(X.'AT) 



E 



3-r(TAT) 



rAT, 



< E 



whence with B" we have found a set that can play the role of B in the 
Lemma's statement. 

Finally, we need to consider the case where T = +00. One has 

e-'-(«An)^(^^^^) = X{«<+oo}e-^('^^")5 (^.An) e'^^C^.) as n ^ 00 P-a.s. 

for all X- valued v that are almost surely finite (for, due to = A 

and g (X+00) by definition of the cemetery A, one has e~'''^g{X.J) = 
X{v<+oo}^~^^ 9{Xv)) ■ By our assumption of uniform integrability, we even 
have L^-convergence in the previous convergence assertion (cf. e.g. Bauer 
[HI Satz 21.4]) and therefore obtain 



lim E 

n^oo 



-rivAn) 



g{X^An) =E[e-™5(^t. 



Combining this result with the Lemma's statement for T < +00 (which has 
been proven before) we get the Lemma's estimate for T = +00, too. 

□ 

Corollary 2.1 (Formula for an option price using hitting times). 

Let X be a d- dimensional basket with an associated risk-neutral measure P 
and discount rate r > 0. Consider a discrete subset Z C [0,+oo). Suppose 
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g = {K — f)y 0, and assume that the process e ^' f{X) is a F-submartingale. 
Then one has 



sup E 

acM'*x[0,T] measurable 

sup E [ 

r Xr\[0,T]-vatued stopping time 



G 

g{Xr 



for all T < +00. If the the set {e ^"^g [Xr] 
uniformly F-integrable, then the equation 



T Z-valued stopping time} is 



sup E 

ScM'^x [0,+oo) measurable 



sup 

r I-valued stopping time 



E[e-'^g{X^)] 



holds. 



Without going into detail we remark that similar results can be obtained for 
non-discrete T as well, as proven by N El Karoui |15| . 

Definition 2.4. Let I C [0, +00) be discrete, G cR'^x [0, +00) and GcR'^ 
measurable, and X a d-dimensional Markov basket with an associated family 
of risk-neutral measures ¥' and discount rate r > 0. We define 



Vlx--{x,t) 



^rt^{x,t) 



ix,t)^g, 
ix,t) £ g. 



as well as 



G,X 



: X 



e-'--GgiXr 



9{x), 



xeG. 



If G I, then we can simply write 



''9{X,. 



and 



VS,x ■■ X 



E^ 



Instead ofVgxix,0), we shall often simply write Vg-^ix). Also, the 
subscript X will be dropped when no ambiguity can arise. Furthermore, Vg 
and Vq will be shorthand for Vg^° and Vq^° , respectively. 

As another notational convention, let us from now on use supg(-]gdx[o,T] ^nd 

SUPGcRrf to denote SUpg(-iRdx[0,T] measurable SUPccRd measurable, respectively. 

The following Theorem is also a classical result from the theory of optimal 
stopping and Snell envelopes (cf. e.g. Neveu (221 Proposition VI-2-8] and El 
Karoui Theoreme 2.31]) 
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Theorem 2.1 (Optimality of the immediate exercise region). Let X 

be a d- dimensional Feller basket with P' being an associated family of risk- 
neutral measures and r > being the discount rate belonging to P' . Suppose 
g = {K — f) y 0, I C [0, +oo) is discrete, and T £ [0, +oo]. Assume, 
moreover, that e~^' f{X ) is a W -submartingale for all x . Define 



-x,r 



(x, t) G M'' X [0, T] : sup V^{x, t) < g{x) 



F 



if T < +00 (we may drop the superscript T wherever this is unambiguous) 
and else 



F 



Then 



'^'+°° = ) (x, t) eR'^ X [0, oo) : sup V^{x, t) < g{x) 

I gcM''x[0,+oo) 



Vx E R'^ y/i,T (x, 0) = U^{T) (x) 



ifT < +00, and V^i,+oo (a^, 0) = {-\-oo){x) for all x € R'^ such that the set 
{e~'^'^g (Xr) : r Z-valued stopping time} is uniformly ¥^ -integrable. 

Proof. Let T < +00. Using Corollary 12.11 and recalling the definition of 
If^^", all we have to show is 



Vx G 



sup e(^'°) 

x[0,T] 



(where we recall that fg < Tg denotes the first nonnegative entry time into 
Q). However, after exploiting the special particular shape of F^ , we can - 
due to the boundedness of g < K which yields Vg < K for all Q which allows 
us to apply Lebesgue's Dominated Convergence Theorem - swap sup and E 
to get for all x G M"^, 



e-"V sup if (X-i 
gcM'*x[0,T] ^ -p' 



e-'"^^ sup e'"^i E (^"^^ '"^^ / ^-rr- 
gcR'*x[0,T] 



sup e(^'°) 

x[0,T] 



(where for notational convenience Tg should denote the first nonnegative 
entry time into Q). Now, let us use the strong Markov property of the Feller 



16 



process X, and for this purpose, let 9 denote the shift operator on the space- 
time path space D ([0, +cxd), M'^ x [0, +00)) (which is the set of all cadlag 
functions from [0, +00) into x [0, +00) - recall that all Feller processes 
have a cadlag modification). We obtain 



V^x(x,0) 



sup e(^'°) 

gcK'*x[0,T] 

sup e(^'°) 

gcK'*x[0,T] 

sup e(^'°) 

gcK'*x[0,T] 



]E(^o,0) 



But e ''■(^(X.) is a P^'-supermartingale for all x G M°', therefore by Doob's 



Optional Stopping Theorem, (e "^^giX^y) 



ve{f^AT',f^oeTf^AT'} 



must also be 

a P^-submartingale for all x € M'^ and T' G (0, 00) (note that fg <fgO dr^j 
a.s. because of the fact that 6 is the shift operator for the space-time 
process {t,Xt)t>o, rather than simply for X). Letting T' tend to infinity, 
we can employ Lebesgue's Dominated Convergence Theorem (atS g < K 
yields e-™c/(X„) < G ^^(P^) for v e {r| A T', f| o ^-i^ A T'} for all 

T' G (0, +00) and x G M'^) in order to get that the expected value of 
of e"''^s ( X-i ) . Hence 



is always greater than or equal to the expectation 



lfx(x,0) 



sup E(^'°) 
gcR''x[0,T] 



e ''^og[X,, 



> sup e(^'°) 

£;cIR'*x[0,T] 

sup e(^'°VJ(x,0). 

ecR'^x[0,T] 

The case T = +00 can be dealt with analogously. 



□ 



Remark 2.1. Whenever Lemma \l/A mav he applied, the immediate exercise 
region of an American option is the intersection of the immediate exercise 
regions of the corresponding approximating Bermudan options: 



{x,t) G X [0,T] : sup vg 



Vi'^^°'\x,t) < gix) 



pj jpsNo,T ^ 1^ ^2 
s>0 mGN 



17 



and its optimality follows from the general theory of optimal stopping (cf. 
Neveu Griff eath and Snell ;i8J). 

Lemma 2.2 (Time-stationarity of immediate exercise regions for 
perpetual Bermudans). Let X be a Levy basket withF' being an associated 
family of probability measures and discount rate r > 0. Then for all s > 
we have 



for all X 



satisfying the condition that 



{e '"^g {Xt) '■ T sH^-valued stopping time} is uniformly ¥^-integrable. 

Proof. Consider an integer n G Nq, and an x G U.'^ such that 
sup {e~^'^ g (Xt) : r stopping time} is P^-integrable. Then we shift the time 
scale by ns to get 



^rns^{x,ns) 



sup e ''^gglXr- 

SCK<*x[0,+oo) ^ 

sup e~'^^ig (Xt-^ 

gcK'*x[ns,+oo) ^ * 



sup e 

g'CM<'x[0,+oo) 



9 [Xr-g°On~ 



sup e ^' g \ Xr^^^ 
g'ClR'*x[0,+oo) " 



where 6 denotes the shift operator on the space (as opposed to space-time) 
path space D ([0, +00), M'^) . Because of the boundedness oi g < K which 
entitles us to apply Lebegue's Dominated Convergence Theorem, we may 
swap sup and E to obtain 



^rns^{x,ns) 



sup e '"^5 (; ( X. 
gcM''x[0,+oo) 



sup e™'E(^'"") 

gcK'*x[0,+oo) 

sup yj^"(x,ns) 



ecK''x[0,+oo) 

for all re G Nq. Thus we conclude 

sup Ff°(x,t) 

x[0,+oo) 



sup Vf°{x,{)) 

x[0,+oo) 



for all t G sNq. If we insert this equality fact into the definition of F**^", 
we see that the condition determining whether a pair (x, t) belongs to F"^'' 
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does not depend on t. On the other hand, by Corollary 12.11 

sup Vf'> (x, 0) = (+00) (x) , 

gcK<*x[0,+oo) 

and the left hand side equals - by our previous observations in this proof - 
the term featuring in the definition of F''^" . 

□ 

Summarising the previous deliberations, we deduce that under the assump- 
tions of the previous Lemmas, the expected payoff of a perpetual Bermudan 
option of mesh size s > equals 



where := := {x G M"' : U'^o{+oo)[x) < g{x)} 

xGM'^ : sup ^cHd W[e-^-hg{Xri^)] <g{x)\. 

measurable ) 

Remark 2.2. As we have already observed in B.emark \2. 1\ under the as- 
sumptions of Lemma MJA and Theorem \2. 1\ the region G* := HneN ^9 ~ 
r\s>o^g (where G^ denotes the immediate exercise region of the perpetual 
Bermudan option with payoff function g and exercise mesh s) will be the op- 
timal exercise region of the American perpetual option with payoff function 
9- 

Furthermore, whenever G G*, then the immediate exercise region for 
the perpetual Bermudan option with payoff function gxG ctnd exercise mesh 
s > 0, denoted G^^g, will again be G (and will thus no longer depend on 

s): This follows, due to G* = flneN " ~ r\s>o^a> /™"^ more gen- 
eral fact that whenever G Gg for some s > 0, then G^^g = G. (To 
prove this fact, simply observe that on the one hand G'^^g 5 Gg, since 
sup ^cRd E- [e-'-^H (xcg) (Xrij)] < sup Hcmt [e"''^^5 (^r^)] < 

measurable measurable 

g{x). Hence, G^^g 5 G, since Gg ^ G by assumption. On the other hand, 
the immediate exercise region G^^g must be inside G, for the payoff function 
Xg9 trivially vanishes outside G.) 

Lemma 2.3. Let us fix a Levy basket with an associated family of risk-neutral 
probability measures ¥' and discount rate r > 0, as well as a region G C 
and a real number s > 0. Then we have 

Vx^G Vi;{x) = e-'-'¥x,^x.*V^{x). 

In particular, using Lemma \2/2\ one has the following equation for the ex- 
pected perpetual Bermudan option payoff: 

Vx G ?7"^°(+oo)(x) = e-^'Fxo-x^ * t/"^°(+oo)(x). 
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Proof. Using the Markov property of X, denoting by 9 the shift operator on 
the path space D ([0, +cxd), M'^) of a Levy process X, and taking into account 
the fact that Tq > (ie Tq > s) in case x G, we obtain: 



Vx G VSix) 



e-'^O^Ssg o (X{r^>s} + X{r^<s} 



□ 



3 Natural scaling for continuity corrections of 
American perpetuals 

Adopting the terminology of Broadie, Glasserman and Kou [10], we shall refer to 
the difference between an American and the corresponding Bermudan options (on 
the same basket and with the same payoff function) as a "continuity correction". 

3.1 Continuity corrections 

Let X be a d-dimensional Markov basket with an associated family of risk-neutral 
measures P' and discount rate r > 0. Consider a (non-dividend paying) perpetual 
Bermudan options with exercise mesh size s > 0, a measurable set G C R'^, 
and a bounded nonnegative measurable function g :M.'^ ^ M>o (the assumption 
of boundedness applies e.g. to the case of put options). The expected payoff 
of this option with respect to the exercise region G for logarithmic start prices 
X G CG is then given by 

Remark 3.1. We shall consider the problem of pricing an American per- 
petual option with payoff function g. Let G he its optimal exercise region. 
Then it is enough to compute the price for an American perpetual option with 
payoff function g-XG- As was shown in Lemma \l.^ if X has a modification 
with continuous paths and g is continuous, then the American price will be 
the limit, as s I 0, of the expected payoff of a Bermudan option with payoff 
function g ■ XG o-nd exercise mesh size, viz. Vq. Since G was the optimal 
exercise region of the American option, G will he, independent of s > 0, the 
optimal exercise region of a Bermudan option with payoff function g ■ XG- 
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We start with the definition of what will be used as a measure for continuity 
corrections. 



Definition 3.1. Let X be a d- dimensional Markov basket with an associated 
family of risk-neutral measures ¥' and discount rate r > 0. Consider a 
bounded continuous payoff function 

g:Rd, ^ ]R>0. We define 



p(s)(x) := limE"^ 

no 



no 

t=s 



whenever this limit exists. 

Remark 3.2. 1. If X has a modification with continuous paths, then 
(readily due to Lemma \l/A and monotone convergence), 



pis) = E- 



sup e '^^G 

.mGN 



sup E' 

m<=N 



hence the limit in the definition of pis) exists, whence pis) is well- 
defined. 

2. pis) provides an estimate for the continuity correction. Since g is 
bounded and nonnegative. 



VxgCG <sup5-E^ 



so 



m 



< supg • pis) on CG. 
G 



(8) 



3. pis), for any s > 0, depends only on G, not on the payoff function g 
itself. 



3.2 Formulae for p(s) in dimension one 

Formulae for continuity corrections for barrier options in the one-dimensional 
Black-Scholes model have already been derived by other authors (e.g. Broadie, 
Glassermann and Kou 0, Horfelt jSS], Howison [23J, as well as Howison and 
Steinberg [22j). We will nevertheless, for the sake of illustration, show how our 
approach applies to this (compared to other examples, of course, very simple) 
setting. This should be seen as a motivation for the proof of Theorem 13.21 and 
its Corollaries, where higher dimensions are studied. 

Theorem 3.1. Consider a 1-dimensional Levy basket X with an associated 
family of risk-neutral measures F' and discount rate r > 0. Let G = (—00,7) 
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for some 7 G M. Then one has for all s > the relations 



n=l 



-rntTi 



n-1 



exp 



E 

n=l 



fl {X^t > 7} n {Xnt < 7} 
V^nt <0} 



1=1 
„—rnt 



tio 



t=s 



(9) 



(and the limit on the right exists whenever p{s){'y) is well-defined). 

Corollary 3.1. Let G = (—00,7] or G = {—00,^) for some 7 G M and as- 
sume {Xt)t>o = (^0 + ■ -^t + ( ^ ~ ^ ) ) J words: X is the logarith- 



mic price process of the one- dimensional Black-Scholes model with constant 



volatility a and discount rate r > 0. Then, whenever fi := r 



> 0, there 



exist constants cq, Co > such that for all sufficiently small s > 0, 

CQS^ < P(s)(7) < Cos^ 



If both fj. < and r > there exist constants ci, Ci > such that for all 
sufficiently small s > 0, 



cisv^ < p(s)(7) < CiS2. 

Remark 3.3. Although computing the constants Co,Co,ci,Ci explicitly is 
possible, we refrain from doing so for the moment, as it is not required to 
find the right scaling for an extrapolation for Vq from s > to s = and 
it would not provide any additional useful information for our extrapolation 
purposes. The same remark applies to all examples and generalisations that 
are studied subsequently. 

Proof of Theorem AS. 11 The Theorem follows from a result by Feller p. 
606, Lemma 3] on processes with stationary and independent increments. 
For, if we define 

n-1 



n {Xis > 0} n {Xns < 0} 



.i=l 



Vs>0 V(ZG[0,1) :=5]g"P° 

n=l 

then Feller's identity ^Hl P- 606, Lemma 3] reads 

00 ^ 

Vs>0 Vge[0,l) -In(l-e(g,s)) = V— P°{X„, <0} (10) 

z — J ri 



n=l 



and holds whenever X has stationary and independent increments, in partic- 
ular for all Levy processes (note that our definition of a Levy process requires 
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them to be Feller processes in addition). This entails 

Tn-l 



oo 

rnsmO 



e 

n=\ 



1=1 

oo 



l-exp P°{^n. <0} , (12) 



n 

n=l 



which is enough to prove the second identity ((HJ in the Theorem. □ 

Proof of CoroUarv \S.ll First, it makes no difference whether we consider 
Gi = (—00,7) or Gi = (—00,7], as Gi\Go = {7} which has capacity zero. 
In light of Theorem we shall show that if /x > 0, there exist constants 
Co, Co > such that for all sufHciently small s > 0, 

1 / g-r-ns \ 

COS 2 < exp - ^ ^T^"^ < 0} < Cos 

\ n=l / 

2 

and if both /U < and r > there exist constants ci, Ci > such that for 
all sufficiently small s > 0, 



cisV2 < exp ( - ^ ^^"^ ^ 0} ) ^ C'is2 . 



n=l 



Now, the scaling invariance of Brownian motion yields for all n G N and 
s > 0: 

= (2^)"^/^ / exp ^ dx. (13) 



We divide the remainder of the proof, which will essentially consist in finding 
estimates for the right hand side of the last equation, into two parts according 
to the sign of /i. 

Case I: fj, >0. In this case we use the estimates 

I 1 9 9 9 

o 9 \x + v\ y X 

\/x<Oyy<0 -y^-x^<- < , 

- y - y - 2 ~ 2 2 

thus 

Vy < 

e-y" f e-'"dx < f exp f-^±^^ dx < e~y"/^ C e'^'/^^x, 
J —00 J —00 
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hence by transformation for all y < 



vr _„2 f x'^\ , /vr 



Due to equation ifT^. this entails for all n G N, s > 0, ;U > (if we insert 
_M(^s)i/2 for y) 



2^/2 

Therefore for arbitrary r, s > 



< IP° {Xns < 0} < 



* n=l n=l n=l 

The sums in equation have got the shape of ^ for g < 1. Now one 
performs a standard elementary computation on this power series: 

oo „ CO ,.q i-q f? 1 

= E/ ^"^^=/ E'''dr= --dr = -ln(l-g)(16) 

n=l " n=0-^0 n=0 ^0 ^ " ^ 

which immediately gives 



(CO — j'YlS 
-^^^^P°{Xns <0} 



< 1 - e 



-.1 r+ii| 



1 

2V2 



when applied to equation (|T5|l . Due to de I'Hospital's rule, the differences 
in the brackets on the left and right hand sides of the last estimate behave 
asymptotically like s when s | 0. This is sufficient to prove the estimate in 
the Corollary for the case of /i > 0. 

2 

Case II: /U < and r > ^ . In that case we employ the estimates 

Vx<0 V-y<0 < < — - H 

^ - 2 2- 2 - A 2 

and proceed analogously to Case I, to obtain 

< I exp ( — — ) dx < y/ne^ . (17) 



2 - J-oo V 2 
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In the special case of y := ^(ns)^/^ < 0, this leads to the estimate in the 
statement of the Corollary via 



2 2 
P 2ct2 6 2(72 

< F° {Xns <0}< 



V2 



and 



„ — ns r+-^^ „ „ — ns r — ^ 



n=l n=l ^ n=l 

Therefore in case /i = the scaling exponent is exactly ^. □ 

The identity © of Theorem l3.1l can be used to derive estimates in the spirit 
of Corollary 13. II in more general situations. We will illustrate this by means of 
the following example: 

Example 3.1 (Merton's jump-diffusion model with positive jumps 
and "moderate" volatility). Suppose the logarithmic price process X is 
governed by an equation of the form 

yt>0 Xt = Xo + at + f3Zt + aBt 

where a £ E., (3,a > 0, Z is the Poisson process (thus, in this setting, 
only positive jumps are allowed for simplicity) for the parameter 1 and B a 
normalised one- dimensional Brownian motion, and the stochastic processes 
B and Z are assumed to he independent. Let P' he an associated family 
of risk-neutral measures and r > the discount rate. In order to employ 
10), we shall compute the sum 

E~=o {Xns < 0} for all s > 0. Since 

P° {Xns < 0} = P'' < o| for arbitrary n, s we may without loss of gen- 
erality take (3=1. Given a < r — lnE*^Zi, the process exp (Xt — rt)^^Q will 
be a martingale for the unique a satisfying 



1 =E° 



^2 



exp ( a - r + — + Zi 



Let us, in addition, assume a > 0. Then exp {Xt — rt)^~^Q being a martingale 
implies 



0- = y 2 (r - a + InEO [e^i]) < ^2 (r + InE^ [e^i]). 
Now, by definition of the Poisson distribution together with the symmetry 
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and scaling invariance of Brownian motion 

°^ p—rns 

V p°{x„, <o} 



n 

n=0 

n=0 fc=0 

= EE'-^-^ {b, < -^in.)i - tin.ri] (19) 

n=0 fc=0 i, 

V 

= (27r)-l/2 J_^^ " exp(^^jdx 

(with the convention that 0^ = 1). Now let us first of all try and find 
estimates for the probability in the last line. By equation |77P applied to 



'.{ns)h+^(ns)-^Y 



< F°<Bi< (ns)2 (ns)"2 



2V2 I ^ a 



e 

< — 



2 

which yields, using the abbreviation a' := ^, 

-a'^ns-2^k 1^ r , 

ns 2 (ns 2 



< F^\Bi< 



2V2 I cr^ ^ CJ 



SO 



e 2 <T 2cr2ns 
< 



< P° <j < --(ns)5 - -(ns)~5 J> (20) 



2^2 I ' a 



;2 / 

< - 



2 

Thus, we can perform the following estimates to derive an upper bound of 
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the sum in 

\k „—rns 



EE--^-^r«{B.<-^(™)i-f(™)-i} 



n=Ofc=0 

/ ,2 

^ ^ ^k^A' ^ -"^ 

n=0 k=0 



CXJ 



-nsi 



2 ^ n 

n=0 



/ ,2 -q/\ 

2 ^ n 2 \ 

n=0 \ / 

where the last line uses that ^ = > a - [r + ¥P [e^^] ) > and we need to 

impose the condition that ^]) < 1 + r + ^ (which, given r > 

and a, will be satisfied if a > is sufficiently small) to employ the identity 



.21) 

n \ — n 

n=0 



T/ie lower bound follows simply from 

fc=o ■ ^ ^ 

(for n = recall that 0^ = I in this paragraph by our earlier convention) as 
this entails (when exploiting the estimate \2(1\) and finally \21]) ): 

n=Ofc=0 ■ ^ ^ 

2\/2^ n 2^/2 V / 

a consequence of these estimates and using the Taylor expansion of exp 
around 0, we now get the existence of two constants C3 > and C3 > 
(which can be computed explicitly) such that for all sufficiently small s, 

1 / g-rns \ 

C3 • S2 < exp - po {Xns < 0} < C73 . 

V n=0 / 
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Finally, we may apply identity 0) from Theorem I.V. 1\ - as this is an imme- 
diate consequence of Feller's identity fl6\ p. 606, Lemma 3] - and conclude 
that if G = (—00,7), then 

C3 • < p{s){n) < Cs • 
for all sufficiently small s > 0. 

3.3 A Wiener-Hopf type result in higher dimensions 

The proof of Theorem ITU relies heavily on Feller's result [El p. 606, Lemma 3] 
which in turn is proven by means of elementary Fourier analysis and a so-called 
"basic identity" [El p. 600, equation (1.9)]. 

Hence, if one aims at generalising Theorem 13.11 to higher dimensions, one 
should first of all find a multi-dimensional analogue of the identity lfT?Hl. 

Indeed, we shall see that this is feasible. Let us for the following fix a stochas- 
tic process X = {Xt)t>o on with stationary and independent increments. 

Lemma 3.1. Suppose H is a measurable subset ofM.'^, and s > 0. Define 
for allnGN 



fl {x^, G Zh} n {Xn, e K n H} 

l<i<n 



yK eB [R'^J Rn{K) := P° 
as well as 

yK eB (m'^) Qn{K) := P° 



fl {Xis e Zh} n {Xns Zh n K} 



l<i<n 



(in particular Rq = 5o[- H H] =0 and Qq = (5o [• H C//] = 60). Then for all 
n£No, 

Qn+l + Rn+l=Qn*^\,- 

Proof. Consider a measurable K CR'^. Clearly, 



(Qn+l + Rn+l) (K) 



fl {X,s G ZH} n {X^n+l)s G K} 



.1=1 



(22) 



On the other hand, since X is a Markov process, we have 

QniK) = iPsiXCH-)rXKiO) 

(where {Pt)t>o '■= {^Xt * ')t>o translation-invariant Markov semigroup 

of transition functions for the process X whose increments are stationary and 
independent), thus 

' f{y)Qn{dy) = {Ps{xcH-)rfio) 
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for all nonnegative measurable functions /. But this implies 



= iPs{XGH-)r (0) 

= (Ps {XGH-)r ° i^'x. * Xk) (0) 
= iPsiXGH-)r°PsXKiO), 

and the right hand side of this equation coincides with the one of identity 

lEai- □ 

Applying Fourier transforms we obtain 

Corollary 3.2. Let us adopt the notation of the preceding Lemma and define 
the Fourier transform of a countable sequence (^n)„ of finite measures on M'^, 

denoted by Jpnln = /I : (0, 1) x M'^ ^ C, by 

oo „ oo 

Then for all q € (0, 1), and € M'^ the equation 

1 - (iU.(g,C) = iOXiQX) (l - q^siO) 

holds. 

Proof. The result of the previous Lemma reads 

when we apply the Fourier transform. After multiplication with q'^'^^ and 
summing up over n G Nq, one arrives at 

VgG(0,l) VCGM^ C) "XtoW +%, C) "XhW = C)P°J (C), 

=Qo(C) =Ro{C) 

hence 

Vg G (0, 1) R{q, •) - 1 = qQ{q, " Qil, O- 

This is our claim. □ 

Definition 3.2. A subset A Q is called +-closed if and only if A is 
measurable and A + A C A, that is sums of elements of A are again elements 
of A. 
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Lemma 3.2. Let H be dosed and convex with G dH. Whenever both H 
and its complement Ci? are +-closed, there exists a yu € such that 



i7 = |x G M'' : ^xyH > o} . 



Proof. According to the projection theorem, there exists a G M'' such 
that 

ii" C |x G M'' : ^xyu > o} =: K. 

We must show H D K. 

First, note that there is no y G ZH such that ^yyn > 0. For, if there was 
one, then ^{—y)yH < 0, hence —y G C ZH, which yields, because Zh IS 
+-closed, also 

= {-y)+yeZH, 
contradicting our assumption G dH C H. Therefore we already have 

ii" D |x G M"* : ^xyH > o} . 

Now (9{xGM'^ : *xyH > O} 5 {x G M'' : *xyH = 0} (even equality 
holds, but this is not needed here), since whenever ^xyu = 0, then 
j^yn + X — >xasn^oo and also * {^yn + x) yn = ^lynyH > for 
all n G N, making x a limit point of |x G M*^ : ^xyn > O}. 

But on the other hand, H is closed by assumption and we have already 
proven HD{xeR'^ : ^xyn > O}. So we get in addition 

D |x G M'^ : *xyH = o} 
and thus H D K as claimed. 

□ 

Lemma 3.3 (a la Feller, Wiener, Hopf). Suppose H is a +-closed set 
and its complement Cif is a +-closed set as well. Assume furthermore ^ H 
(ensuring Rq = 0), and let In be the main branch of the logarithm on C. Then 

- In (l - R{q, O) = E - / (dx) 

for all {q, Q G (0, 1) x such that the left-hand side is well-defined. In 
general, for all q G (0, 1), one has at least 

l-i?(g,0)=exp -J] ^ ^ 



n=l 



30 



Proof. Let q £ (0, 1). According to the previous Corollary 18.21 we have 

VCGf/ In i =ln i lnQ(^,C) (23) 

wherever this is defined. Due to the identities Yl'^=i IT ~ 

r £ Bi{0) C C (cf. equation (fTH|l in the proof of Theorem 13.11 above) and 

P'^Xs = IP'^Xs*" this can also be written as 



e^-*C^(P%J*^rfx) (24) 

n=l -^IK" 

f:^(«(''.o)%i:ti^(Q(.o-i)". 



n \ / — ' n 

n=l n=l 



However, at least for C = and arbitrary choice of q, one may still state 
identity (|23|l as this follows from Corollary 13.21 more or less directly: First 
we note that 



ln(^(^l-gPOx.(g,0)j •(Q.)J(?,0) 
In (oXil^ 0) - In f 1 - q^^ {q, 0) 



(as in these statements the arguments of In are positive, hence surely in the 
domain of In) and written in series notation 

OO oo ^ 

n=l n=l 

= E^((i-9^(^'0))-(^(9,o)-i 

n=l 

But Corollarv 13.21 implies 

f_l)n // — \ \n 1—— — 

VnGN l-gpo^^(g^o) •(Q„)„(g,0)-1 =-(i?„)„(g,0). 



n w / / n 

Combining these two equations yields (|23|) . Next, note that 

oo 
n=0 

is still a finite measure - concentrated on if - and thus possesses a Fourier 
transform. Analogously, the measure nq^g := Yl'?^=o1^Qn concentrated on 
CH and also has a Fourier transform as it is finite. Now, for arbitrary n G N, 
the properties of the Fourier transform imply 



1 ' 

{Q{Q, •) - l) = [f^Q,q - ^o) = ((/iQ,g - ^( 
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But since H and CH are +-closed sets, \e H + H (1 H and Cff + C-ff C 
Ci? , the measures on the right hand sides of these two equations, ^J'R,q*^ 
and (/iQ,q — Sq)*^, have to be (signed) measures on H and CH, respectively. 
Let us now spht the sum in H24p and insert the terms we have previously 
identified: 

VC G [/ 

n=l n=l -^^^ 

= E-W")(0 + E^(^Q.«-'^o)*"(C). (25) 

n=l n=l 

It is the injectivity of the Fourier transform that yields from this 



_ n \ / — • n 

n=l n=l 

= E ^ (/^^./'^) + E 4^ (/^Q.. - ^o)*" • 

n=l 71=1 

Either side of this equation equals the sum of two (signed measures), and we 
recall that the first measure on the left hand side and first measure on the 
right hand side are both concentrated on H, whilst the second measure on 
the left hand side as well as the second measure on the right hand side are 
both concentrated on ZH. The only way for this to be true is that the two 
measures that are concentrated on each of H or Ci/ are equal: 



oo ^ oo . 

E!-((rVr(-niJ)) = E^(™ 



and also 



n \ / ^ — ' n 

n=l n=l 



E {r-S'' (• n CF)) = ^ ^IL (^^^^ _ s,r , 

71=1 n=l 

the former identity being exactly what the statement of the Lemma expresses 
in the language of Fourier transforms. □ 



■,q ) 



3.4 Continuity corrections in higher dimensions 

On the basis of Lemma IT3l we may now partially generalise Theorem 13.11 to 
higher dimensions when we require G (the set that we refer to the exercise 
region) to be +-closed set. 
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Theorem 3.2. Let X be a d- dimensional Levy basket with an associated 
family of risk-neutral measures ¥' and discount rate r > 0. Also, consider a 
G of the shape G = 7 + if for some 7 G M*^ and some measurable set H such 
that ^ H, and such that both H as well as C-ff are -\--closed. Then for all 
s > 0, 



n=l 
00 

E' 

n=l 

exp 



-rni Ti 



-rntmfi 



n-l 



fl G Cg} n {Xnt G G] 

i=l 

fl {Kit G n {Xnt G H] 



n=l 



4 = 1 



{Xns G G - 7} 



(and the limit on the right exists whenever p{s){'j) is well-defined). 

Proof. The Theorem follows directly from Lemma, 13.31 — in the same man- 
ner in which Theorem 13.11 followed from Feller's original result (ini p- 606, 
Lemma 3]: For, the second equation in Lemma IT3l mav be read 



n=0° 



exp 



00 n 

q 



E 

n=l 



my 



n 



that is 



n=0° 



fl {x,seZH}r\{XnseH} 



l<i<n 




exp 



for all q G (0, 1), in particular for = e 



□ 



Corollary 3.3. Let X be a d-dimensional Black-Scholes model with discount 



t for 



rate r > and volatility vector a (i.e. {Xi)^ = (Xj)Q + crji?j*^ + 

allt > and independent normalised Brownian motions B^^\ i E {1, . . . ,d}). 
Let ¥' denote the corresponding family of risk-neutral measures. Also, con- 



sider a region G of the shape G = 7 + G 



ax 



<0} /< 



or some 



2 ' 



£jL. 
2 



a,7 G M . Let 5 := y X]j=i l^i'^il 0'''^^ A* •= 

Then, whenever ^a^ > 0, there exist constants cq,Co > such that for 
all sufficiently small s > 0, 

CQS^ < /9(s)(7) < Cqs^ 
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I'o 1^ 

// both < and r > > ^^^re exist constants ci, Ci > such that for 

all sufficiently small s > 0, 

1 1 

ClSV5 < /9(s)(7) < ClS2. 

Proof of Corollaru \HJ-A Observe that 

But since the components of X are independent multiples of Brownian mo- 
tions with linear drift at rates /ii, . . . , /i^, respectively, the process 
is a multiple of a normalised Brownian motion with linear drift at rate *a/x: 
For each t > 0, the random variable 

is distributed according to 

"^Lxt = ^^a^,■t,Y.tl\'^^<'^?t = ^*a/..t,<52f (26) 
Hence for any normalised 1-dimensional Brownian motion S, one has 

and 

We can now apply the estimates from the proof of Corollary 13.11 to the 
process {pBt + ^afit)^^^ en lieu of what is X there. 

For instance, if > 0, we shall employ the bounds 



U H= it 

e e ^ 
<F"'{XnseG}< - 



2^/2 

and if *a/i < 0, we will use the estimates 



It |2 It |2 

— - — TTT—ns - — -J—ns 

e 2^2 g 2«2 



□ 
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Appendix 



Reformulation of the perpetual Bermudan pricing 
problem in and 

A Non-applicability of the 
L2(M^) Spectral Theorem 



Consider a d-dimensional Levy basket X with associated family of risk-neutral 
probability measures P' and discount rate r > 0. 
Fixing h > and defining 



we can rewrite the result of Lemma lZSl as follows: 

PA (V^ - g,) = -PAgi (27) 



where we assume that g has a square-integrable extension from G to the whole 
of W^; given this assumption, the gi G L?{W^) of the previous identity can be 
any such extension. 

We will suppress the superscript of A for the rest of this paragraph. 

Also, without loss of generality, we will assume in this Chapter that the 
components of the basket X when following the Black-Scholes model all have 
volatility 1. 

Lemma A.l. Let X he a Levy basket with associated family of risk-neutral 
probability measures V and discount rate r > 0. Then A and PA \ L^(M'^\G) 

are invertible. Furthermore, the L^ norm of A is bounded by (l + e~^^^ ^ , if 

X = n - +B (thus // = (r — ^)'^_-^) where B is a standard Brownian motion. 
Moreover, A is a contraction if fi = 0. 

Proof Suppose u e L^{R'^\G) and u is bounded. Then ess sup \u\ / 
and we may choose a set -ff C G of positive Lebesgue measure such that 
e~'''*ess sup |n| < for all x ^ H (this is possible because r,h > and 

therefore e~'''^ < 1), we deduce 



yxeH 



< e ^^ess sup \u\ 



< u{x), 

which means that 

\/xeH PAu{x) = u{x) - e-''^Fxo-x^ * u{x) ^ 0, 
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hence PAu / (for H has positive Lebesgue measure). So 

ker PA \ L^{R'^ \G) = ker PA n ^^(M^ \G) = {0} 

and we are done for the invertibihty of PA \ ^^(M'^ \ G). Similarly, one 
can prove the invertibihty of A. Finally, A is seen to be a contraction by 
application of the Fourier transform: The Fourier transform is an isometry 
(by Plancherel's Theorem), thus 



(<5o-e-'-¥xo-x,)*/ 
(<^o-e-'-Vxo-x,)*/) 



L2(Rd) 



L2 



-rh+ih^ fix- 



f 





2 




2 \ ^ 


2 






dx 



Now, the factor in front of 

,2 



f{x) in the last line can be bounded by 

(l + e~^'*)^, and it is strictly less than one for /x = 0. Using Plancherel's 
Theorem again, this yields the result. □ 



Now, this is sufficient to apply a Wiener-Hopf factorisation (for a general 
treatment of this kind of factorisations, one may consult e.g. Speck (40], our 
application uses in particular 03 11. Theorem 1]) and state 

Theorem A.l. Let G C R'^ and let X he a Levy basket with associated 
family of risk-neutral probability measures F' and discount rate r > 0. Then 
Vq, the expected payoff of a perpetual Bermudan option for G with exercise 
mesh size h > and payoff function g, is - using the above notation - given 
by 



VA 



\G)) 'pAgi 



gi - A-^PAZ'PAgi 



= 51 - [PA \ L' 

where A = A^A^ is a Wiener-Hopf factorisation of A. 

We observe the following: 

Lemma A. 2. The Hilbert space operator A : L^{R'^,C) 
normal. 



-1 



L2(M^, 



is 
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Proof. We define p := e — ^dxd^'^ (where A*^ is the d- dimensional Lebesgue 
measure) and via the Pubini Theorem one has for every /, 5 G L'^{U.'^) 

{Af,g) = {5o-p){x-y)f{y)dyg{x)dx 

= fiy) i^o -p){x- y)g{x)dxdy 

= (/, {6-po{-I)*g)), 

that is 

A* = {So-po{-I))*- 
But since the convolution is associative and commutative, this implies 

A*A = (5o-po(-I))*(5o-p)*- 
= {So-p)*{6o-po{-I))*- 
= AA*. 

□ 

However, it will not be possible to find a basic system of eigenvectors and 
eigenvalues for this operator, since 

Lemma A. 3. The operator A fails to be compact. 

Proof. Any normalised basis provides a counterexample for the compactness 
assertion. □ 

Therefore, the equation (j27|l cannot easily be applied to connpute the ex- 
pected option payoff by means of a spectral analysis. Thus, our examination 
of the Hilbert space approach in the second part of this Chapter has led to a 
negative outcome. 

However, one can also conceive of the operators A^ as operators on the 
Banach space {W^): 

B The operator equation: 

analyticity in the exercise mesh size 

From now on, /i will no longer be fixed and we will therefore write A^ instead 
of .4. 

If we now assume gi to be an integrable extension of G to the complement 
of CG as an element of Quite similarly to lA.ll we can prove 
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Theorem B.l. Let G and let X he a Levy basket with associated 

family of risk-neutral probability measures ¥' and discount rate r > 0. Then 
Vq, the expected payoff of a perpetual Bermudan option for G with exercise 
mesh size h > and payoff function g, is - using the above notation - given 
by 

Vh = 9i- [pa \ L\R^ \ G)y' PAg, = 51 - A^'PAz'PAgi, 

where A = A^A^ is a Wiener-Hopf factorisation of A. 

It suffices to observe that A^ is - due to the norm estimate for the 
convolution of two integrable functions (as the product of the norms of the 
convolved functions) - also a bounded operator on L^ (M°'). 

We shall now identify g and gi. 

Theorem B.2. With the notation previously introduced, we define E to be 
the semigroup 

iEt)t>o = {e-'-\t,t * •)i>o = * •)t>o ^ ^ (i^^M'^), ^^(M^)) , 

where 

is the distribution of the logarithmic price vector at time t. Suppose x ^ 
G and, with the notation from the previous chapters, X is a (normalised) 
Brownian motion with (possibly zero) drift fBlack-Scholes modelj. Then 
h ^ Vq{x) is real analytic in h on (0, +co) as function with range in the 
Banach space L^iW^). 

Proof. It is obvious that is a semigroup. According to [El Theorem 1.48], 

the set 

£: := | / g L^{M.'^) : t ^ Etf £ L^{R'^) entire} 

is dense in L^(M'^). Hence it is possible to approximate every 5 by a sequence 
{9k}k C£ in L^{R'^). Since 

Vt>0 < e-^* < 1, 

we obtain Etgk Etg for A; —> 00 uniformly in t on M+, where Etgk is entire 
for every t > and k £ N. Thus, t ^ Etg, and thereby t 1— > Pg — PEtg, is 
an analytic function on (0, +cxd) taking values in the Banach space L^iW^). 
Now observe that for arbitrary open U CC M_|_ (the symbol "CC" indicating 
that U is contained in a compact subset of M>o) the following equations 
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hold: 

Vt G [/ = g-(l-PEt\ PL^W") ) \Pg- PEtg) 



-1 



= g-[Y,{PEtf\{Pg-PEtg) 

\fc=0 / 

oo oo 

= g-Y,{PEtfPg + Y,{PEtf^^g 

k=0 k=0 

oo 

= Y^iPEt)' (g - Pg) , (28) 

fc=0 

since the sums converge uniformly in t on t/ CC M4-, yielding the analyticity 
oi t Vq as a function whose range lies in the Banach space L^(M'^). □ 

Lemma B.l. Let u > 0, n £ N. Then the equation 



Eu= [-r + -A- VV gt^,,t * ■ 



holds (where *y denotes the transpose of a vector y). In particular, if ^A/ — 
AtV/ = A/ for some A, /, 

— i^„/ = (-r + Ar5,,,*/. 
Proof. According to Davies jT^ Proof of Theorem 2.39], we have 

An 

—E^={ZE^/^)\ (29) 

where Z denotes the infinitesimal generator of the semigroup E. Now, define 
C to be the convolution operator semigroup {gtfj.,t * ■)t>o (normalised) 
Brownian motion with drift n (as before denoting by 5^^0-2 the Lebesgue 
density of the Gaussian distribution centered around z of variance o"^ for 
all z G M"^ and o- > 0). It is well-known (cf. e.g. [33 p. 352]) that the 
infinitesimal generator of this semigroup C is 

L := ^A + VV. 

By our requirements on /, Lf = on [/. Furthermore, L and C commute: 

Vt > CtL = LCf 

Thus, 

Vt>0 ZEt = -Et = ^ie-'' -Ct) 
d dt^ ' 

= -re-^'Ct + e-^'^Ct 
dt 

= e-^'Ct{-r + L), 
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which due to equation l(29|l already suffices for the proof of the Lemma in 
the general case. And if / is an eigenfunction of L for the eigenvalue A, one 
has (-r + L)"/ = (-r + A)"/. 

□ 

Theorem B.3. The Taylor series for the expected payoff of a perpetual 
Bermudan option as a function of the exercise mesh with respect to a fixed 
exercise region G is for all s > 0: 

oo oo / m ^ \ 

(?!,..., U G No'" 
(^XRd\G ■ * •) (^-r + + ^ ixcg) , 

where, in order to avoid confusion with pointwise exponentiation, denotes 
for any operator A. 

Proof. We know about the real analyticity of t i—^ Et on M.^^ and even, 
thanks to the previous Lemma, the explicit Taylor series. Thereby we also 
have the Taylor series for 1 1-^ PEf. So we can use equation l(28|l and see by 
means of a binomial expansion 

CO oo / oo / \ ''^ 

k=0 k=0 \ i=0 I 

oo oo ™ \l 

n=Om=l + + = „ i=l 

oo oo 
n=0 m=l 

m ^ 

h-\ \-lm = n 



m 



□ 



This Taylor series fails to provide any straightforward possibility for the com- 
putation of Vg- Instead we state the following immediate Corollary of equation 
(f2R|l : 
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Corollary B.l. With the notation as in the previous Theorem, 
, ,00 

V5>0 -VS = ^E«'''"^'(xm<^\g-(5w*-))°"(xg-5)- 

m=l 
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